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ABSTRACT 


Re-entry  of  a  rotating  symmetrical  missile,  which  is  assumed 
to  move  along  a  straight  path,  is  examined.  The  equations  of  rotational 
motion  are  reduced  to  one  second  order  differential  equation  for  the 
angle -of- attack,  the  other  angles  being  then  obtainable  as  quadratures. 

This  form  of  the  equations  of  motion  is  suitable  for  numerical  integration, 
as  all  of  the  exact  constants  of  motion  have  already  been  integrated  out. 

Small  angle  oscillations  are  considered,  and  it  is  shown  that  pre¬ 
vious  analyses  of  the  effect  of  rotation  on  oscillation  are  in  error,  due 
to  an  improper  procedure  for  obtaining  the  small-angle  equations.  It 
is  pointed  out  that  the  precession  rate  of  the  oscillational  motion  is  likely 
to  depend  as  much  on  non-linearities  in  the  aerodynamic  restoring  torque 
as  on  the  rotational  velocities. 

For  a  rotating  missile  with  large  initial  angle-of -attack,  the  adia^ 
batic  invariant  is  used  to  calculate  the  amplitude  of  oscillation  at  the  lower 
altitudes. 
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RE-ENTRY  OF  ROTATING  MISSILES 


by 

Conrad  L.  Longmire 


I.  INTRODUCTION 


In  this  note,  we  discuss  the  angular  motion  of  missiles  which  are 
rotating  as  they  re-enter  the  atmosphere,  without  necessarily  malting  the 
assumption  that  the  angle- of- attack  is  small.  We  do  assume  that  the  tra¬ 
jectory  of  the  missile  is  a  straight  line,  that  the  missile's  external  enve¬ 
lope  is  axially  symmetric  about  one  of  the  principal  axes  of  inertia,  and 
that  the  center  of  mass  lies  on  this  axis. 

II.  THE  EQUATIONS  OF  MOTION 

We  use  the  coordinate  systems  shown  in  Fig.  1.  The  origin  of  the 
coordinates  is  at  the  center  of  gravity  of  the  missile.  The  missile  has 
moment  of  inertia  A  about  the  X'  and  Y'  axes  aind  moment  of  inertia  C 
about  the  Z'  aixis. 

When  the  angle  0  does  not  vanish,  the  aerodynamic  forces  on  the 
missile  lead  to  a  torque  whose  magnitude  T(0  ,  t)  depends  on  0  (and  on 
the  time)  and  whose  direction  is  parallel  to  the  line  ab  in  Fig.  1,  i.  e.  , 
the  line  of  intersection  of  the  X' Y'  plane  with  the  XY  plane.  Thus,  the 
components  of  the  torque  about  the  Z  and  Z'  axes  vanish, 


Tz  =0’ 

(1) 

Tz'  =  0' 

(2) 

and  the  components  about  the  X'  and  Y'  axes  are 

T^,  =  -  T(0  ,  t)  cos  (t>. 

(3) 

Ty,  =  T(0,t)sin0. 

(4) 

The  choice  of  sign  is  such  that  positive  T(0  )  implies 
tries  to  reduce  0  . 

a  torque  which 

Fig.  1.  Coordinate  Systems.  The  XYZ  axes 
direction  of  travel  of  the  missile  in  the  minus 
axes  are  fixed  in  the  missile,  and  the  Z'  axis 
angles  Q  ,  ip  ,  and  <f)  are  the  Euler  angles. 


are  fixed  in  space,  with  the 
Z  direction.  The  X'Y'Z' 
is  the  symmetry  axis.  The 


I 
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Let  p,  qt  and  r  be  the  angular  velocities  of  the  missile  about  the 
X' (  Y'  and  Z'  axes,  respectively.  Then  Euler' s  Dynamical  Equations 


-  (A  -  C)  qr  =  -  T(e  ,  t)  cos  <f> 


A-^-  +  (A  -  C)  rp  =  T(0  ,  t)  sin  ^ 


■cR— 


From  Eq.  (7),  the  angular  velocity  aoout  the  symmetry  axis  is  constant, 


r  =  constant  =  v.  (8) 

We  shall  call  v  the  spin  rate. 

Because  the  torque  T(6  ,t)  depends  explicitly  on  the  time  (on  account 
of  the  changing  air  density),  the  total  energy  is  not  a  constant  of  motion. 
However,  multiplying  Eq.  (5)  by  p,  Eq.  (6)  by  q,  and  adding,  gives 

d  ,2^2,  2  T(e  ,t)  r  A  ■  x  ^  /o. 

(p  +  q  )  - - ^  [  P  cos  (^  -  q  sin  0  ]  .  (9) 

At  this  point  it  is  convenient  to  introduce  the  Euler  angles  instead  of  p, 
q,  and  r.  The  inter-relations  are  (dot  means  time  derivative): 


p  =  sin  0  sin  ^  +  0  cos  (j) 

.  • 

q  =  Ip  sin  0  cos  (ft  -  Q  sin  <f> 

p  =  r  =  (f)  +  Ip  cos  Q  . 


From  Eq.  (12)  one  can  obtain  the  angle  0  once  ip  auid  0  are  known;  but 
we  are  not  interested  in  0  since  it  is  essentially  not  observable.  How¬ 
ever,  using  Eq.  '  8  (10)  and  (11)  in  Eq.  (9)  gives 


[  sin  0  0  +6  ] 


2  T(e,t)  ^ 


1 


If  T(6  ,t)  did  not  depend  explicity  on  the  time,  this  equation  would 
immediately  yield  the  total  energy  constant.  The  left  hand  side  is  the 
rate  of  change  of  the  kinetic  energy  of  rotation,  less  the  spin  energy, 
and  the  right  hand  side  is  the  rate  at  which  the  aerodynamic  torque  does 
work  (a  factor  two  and  the  moment  of  inertia  A  have  been  divided  out). 

£q.  (13)  contains  two  dependent  variables,  6  and  ij).  The  latter  can 
be  eliminated  by  use  of  another  constant  of  motion,  namely,  the  component 
of  angular  momentum  about  the  Z  axis.  That  this  component  is  constant 
is  a  result  of  Eq.  (1).  From  Fig.  1,  it  is  seen  that  the  Z-component 
of  the  angular  momentum,  expressed  in  terms  of  p,  q,  and  r,  is 

=  A  sin  0  [p  sin  0  +  q  cos  (^]  +  C  r  cos  0  (14) 

Replacing  p,  q,  and  r  by  means  of  Eq.  's  (10)  -  (12),  we  find 

2 

sin  0  ^  +  b  p  cos  0  =  (3  (15) 


where  we  have  defined  two  new  constants 

K  -  C 

(16) 

^Z 

P  =  -X- 

(17) 

Eq.  (15)  can  be  used  to  eliminate  ij)  from  Eq.  (13),  with  the  result 

,2 


-at  ® 


2T(0,t)  A  d  ((3  -  bp  cos  0) 

— A~  ~sr  — -7-z - 

sin  0 


(18) 


This  form  of  the  equation  of  motion  in  the  0  variable  will  be  useful  for 
discussing  the  adiabatic  invariant.  For  the  purpose  of  actually  solving 
the  equation  it  is  more  convenient  to  carry  out  the  time  differentiations 
indicated.  The  equation  then  takes  the  form 


t  • 

0 


T(0  ,t) 

“S — 


a 

“STF 


(p  -  b  p  cos  0  )^ 

■"  " .  ir- - 

2  sin  0 


(19) 


This  equation  has  the  form  of  Newton' s  law  for  a  particle  moving  in  one 
dimension  under  a  force  given  by  the  right  hand  side.  The  numerical 
solution  of  such  equations  can  be  easily  accomplished  by  well  known 
methods. 
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Once  0  (t)  has  been  found  from  Eq.  (19)*  ^(t)  is  found  by  quadrature 
from  Eq.  (15). 

There  is  very  little  difference  between  the  analysis  above  and  the 
usual  theory  of  the  symmetric  top  (see,  e.  g.  ,  Osgood,  Mechanics,  sec. 
18).  The  analysis  was  presented  here  to  make  it  clear  that  the  time 
dependence  of  the  torque  does  not  prevent  the  reduction  of  the  equations 
of  motion  to  one  second  order  differential  equation  in  one  unknown. 


III.  SMALL- AMPLITUDE  MOTION  WITH  TIME- INDEPENDENT  TORQUE 


In  this  section  we  shall  assume  that  the  torque  does  not  depend  ex¬ 
plicitly  on  the  time,  and  that  the  torque  iF*ao  large  that  it  dominates  the 
motion  to  such  an  extent  that  0  moves  in  a  small  interval  not  far  from 
0=0.  We  approximate  the  torque  by  a  linear  form,  letting 

■T’  =  “a  ® 

The  meaning  of  the  constant  u  can  be  seen  by  inserting  this  form  in 
Eq.  (19):  to  is  the  frequency  of  the  "pendulum”  type  of  motion  that 
occurs  when  the  angular  momentum  terms  on  the  right  in  Eq.  (19)  are 
not  present.  Thus,  cOg^  may  be  called  the  frequency  of  small  oscillations 
due  to  aerodynamic  forces  alone,  or  simply  the  aerodynamic  frequency. 
tOg.  depends  on  the  air  density,  the  missile  velocity,  the  moment  of  inertia 
A,  and  various  aerodynamic  coefficients  relating  to  the  missile  shape. 

In  this  section,  we  regard  as  a  constant. 

We  now  wish  to  find  the  effect  of  the  angular  momentum  on  the 
frequency.  To  this  end,  we  expand  the  second  term  on  the  right  in  Eq. 
(19)  in  a  power  series  in  0  ,  keeping  terms  up  through  the  first  power  of 
0.  Eq.  (19)  then  becomes 


e  =  -  +13  (P  -  b»/)^  + 


p  b  v]  0  + 


(P-^ 


The  solutions  of  this  equation  are  well  known.  It  is,  in  fact,  the  equation 
of  motion  for  the  radial  coordinate  of  a  two-dimensional  harmonic  oscil¬ 
lator,  with  the  frequencies  equal  in  the  two  directions  and  given  by 


+  (P  -  bv)' 


euid  with  angular  momentum  (P  -  hp)  about  the  center.  If  ^  and  t?  are 
the  (fictitious)  Cartesian  coordinates  of  this  oscillator«  so  that 


the  motion  in  r)  makes  an  ellipse,  as  indicated  in  Fig.  2.  The  fre¬ 
quency  CO  above  is  the  frequency  of  traversing  the  ellipse.  The  frequency 
of  the  ©  motion 


Fig.  2.  Fictitious  Two-Dimensional 
Oscillator 


is  twice  CO,  but  in  order  to  use  a  definition  of  frequency  which  goes  over 
to  the  commonly  accepted  definition  when  there  is  no  angular  momentum, 
we  shall  call  co  the  frequency,  rather  than  2co. 

It  must  not  be  concluded  that  Fig.  2  gives  a  head-on  view  of  the 
trajectory  of  the  missile  nose.  It  is  the  Euler  angle  ip,  rather  than  the 
fictitious  angle  tan“^  Tj/^  ,  which  determines  the  angular  position  of 
the  missile  nose  when  projected  into  the  X,  Y  plane  in  Fig.  1.  The 
aungle  ^  is  to  be  found  from  Eq.  (15), 


I 


^  -  hp  cos  0 
- - 

sin  © 


dt 


(23) 


Into  the  integrand  we  have  to  insert  0  as  a  function  of  t.  From  the 
analogy  of  Fig.  2,  it  is  easy  to  see  that  ©  has  the  form 


e 


(24) 


I 


i 

I 


I 

I 

i 


=  yx  +  sin  2  0)  t 


apart  from  an  arbitrary  additive  phase  in  the  sine.  Here  X  and  /i  are 
constants.  One  relation  between  X  and  jx  is  imposed  by  the  equation  of 
motion,  Eq.  (21);  substituting  Eq.  (24)  therein,  one  can  find  that 


v2  2  v\2/2 

X  -  =((3-bv)  /o) 


(25) 


This  equation  is  equivalent  to  the  statement  that  the  potential  energy, 
corresponding  to  the  force  term  in  Eq.  21,  must  have  the  same  value 
at  the  two  turning  points  0  .  and  0  (see  Fig.  3).  There  is  no 

other  relation  between  X  and  p,  as  the  total  energy  of  the  0  motion  is 
arbitrary. 

When  the  expression  (24)  is  substituted  for  0  in  Eq.  (23),  it  is 
seen  that  the  integral  cannot  be  performed  exactly.  However,  we  may 


V(0) 


Fig.  3.  Energy  Diagram 
for  0  motion. 


expand  the  integrand  in  Eq.  (23)  as  a  power  series  in  0,  keeping  the 
same  number  of  terms  (three)  as  we  kept  in  deriving  Eq.  (21).  Eq.  (23) 
then  becomes 

0 

The  expression  (24)  is  to  be  inserted  for  0. 

Let  us  compute  the  increase  in  ip  during  one  period  of  the  oscil¬ 
lation,  i.  e.  ,  carry  the  integration  from  t  =  -  tt/u  to  t  =  +  tt/w.  One 
then  finds  (with  the  help  of  Eq.  25)  that  the  term  (P  -  b  v)/0^  integrates 
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to  Ztt.  In  fact,  this  term  represents  just  the  angular  rate  that  one  would 
have  in  the  fictitious  two-dimensional  harmonic  oscillator.  Thus,  for 
this  term  alone,  Fig.  2  would  give  the  correct  head-on  view.  However, 
the  other  terms  in  the  integrand  of  Eq.  (26)  lead  to  an  average  rate  of 
advance  of  the  angle  of  aphelion  of 


At  / 
ap 


-  P 
"  T 


hp 


7  bv  . 

-nir)  ® 


(27) 


where  0  ^  =  X  is  the  time  average  of  0 


Because  of  the  expansions  used  above,  all  the  results  of  this  section 
are  valid  only  when 


0 


«  1 


(28) 


By  calculating  the  angle  0  at  which  the  force  term  in  Eq.  (21),  i.  e.  , 
the  right  hand  side,  vanishes,  one  sees  that  the  small  angle  requirement 
demands  generally  that 


CJ  ^ 

a 

2  k2  2 

w  »  b  V 
a 


(29) 


Thus,  generally  the  results  of  this  section  are  valid  only  when  the  aero¬ 
dynamic  frequency  is  much  greater  than  the  angular  velocities  of  rotation. 


The  only  exception  to  this  rule  is  the  case  when  p  is  very  nearly 
equal  to_  h  p  ,  in  which  case  we  do  not  need  large  in  order  to  have 
small  0  ^  A  spinning  object  with  little  wobble  (i.  e.  ,  with  most  of  its 
angular  momentum  about  the  symmetry  axis),  and  with  a  small  angle  0 
between  the  symmetry  axis  and  the  Z  axis,  has  p  very  nearly  equal  to 
b  V.  Thus,  this  is  the  case  of  a  gyroscope  (i.  e.  ,  no  torque)  nearly  aligned 
and  only  slightly  wobbling.  For  this  case,  Eq.  (22)  gives  the  frequency 
u  =  hp/2  and  Eq.  (27)  gives  the  precession  rate  =  hp/2.  Both  of  these 
results  are  correct  for  the  slightly  wobbling  gyroscope. 


The  corrections  to  the  frequency  due  to  rotation,  and  the  precession 
of  the  aphelion,  differ  from  expressions  given  for  these  effects  in  the 
literature.  ^>2  The  reason  for  this  is  that  the  small  angle  approxima¬ 
tion,  in  the  references  cited,  was  not  made  to  sufficient  accuracy  to  get 
the  frequency  correction  and  the  precession  rate  correctly.  In  these 
references,  the  approximation  made  was  to  replace  cos  0  by  unity  and 
sin  0  by  0  .  However,  in  arriving  at  the  expansion  in  Eq.  (21),  it  was 
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necessary  to  take  the  first  three  terms  in  the  expansion  of  the  sines  and 
cosines  (in  the  term  arising  from  rotation)  in  order  to  get  the  first  order 
correction  to  the  frequency. 


Let  us  estimate  the  error  induced  by  neglecting  higher  powers  of  0 
on  the  right  hand  side  of  Eq.  (21).  The  rotational  term  in  Eq.  (19)  (i.  e.  , 

i  3 


the  second  term  on  the  right)  will  contribute  a  term  in  0 

n2 


The  coefficient 

of  this  term  will  be  a  linear  combination  of  b^  v^,  and  pbv.  The  tor¬ 
que  function  T(0  )  is  likely  to  be  non-linear.  If  we  assume  that  T(0 )  is 
expandable  as  a  power  series  and  is  an  odd  function  of  0  ,  the  next  term 
in  T(0  )  would  also  be  of  order  0  Because  the  torque  function  is  likely 
to  be  quite  non-linear,  the  coefficient  of  the  0  ^  term  is  likely  to  be  as 
big  as  the  coefficient  of  the  linear  torque  term,  namely  ,  Thus,  we 
should  consider  the  effect  of  adding  a  term 

3 

K  0 


to  the  right  hand  side  of  Eq.  (21),  where 
K  =  K 


rot 


=  O  (P^.  b^  pbi^) 


for  the  rotational  corrections,  and 
K  =  =  O  («^) 

for  the  non-linear  torque  corrections.  The  latter  correction  is  apt  to  be 
much  the  larger,  since 


0 


w 


according  to  Eq.  (21),  and  we  require  0  to  be  small.  Now,  adding  a 
term  K  0  ^  to  the  right  haind  side  of  Eq.  (21)  will  perturb  its  solution 
by  O(K0  );  that  is,  if  ©  o(f)  and  ©^(t) 

perturbed  solutions,  respectively. 


are  the  unperturbed  and 


ej(t) 


=  0Q(t)  [  1  +  o  ( 


K  es 


U)  ' 


)  ] 


Thus,  the  fractional  correction  to  ©  o(t)  is  of  order  ©q  if  K  =  K  but 
of  order  0  if  K  =  K-p.  It  can  then  be  seen  from  Eq.  (26)  that  file  non¬ 
linear  torque  correction  gives  a  correction  to  the  rate  of  advance  of 
aphelion  which  is  of  the  same  order  as  the  rate  previously  calculated, 
Eq.  (27).  Thus,  a  non-linear  torqne  can  completely  modify  the  pre¬ 
cession  of  the  aphelion  angle. 
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IV.  SMALL- AMPLITUDE  MOTION  FOR  EXPONENTIALLY  VARYING 


We  again  assume  that  the  motion  is  confined  to  small  angles,  but  let 
the  aerodynamic  frequency  vary  with  time  according  to  the  approximation 


(30) 


We  shall  call  the  constant  a  the  "density  rate".  It  is  related  to  the 
missile  velocity,  the  re-entry  angle  and  the  atmospheric  scale  height.  A 
typical  value  for  a  is  about  l/4  sec"^. 

Equation  (21)  covers  this  case,  since  was  not  differentiated  in 

deriving  it.  Let  us  re-write  Eq.  (21)  in  the  form 


•  • 

0 


at 


e 


where  we  have  defined 


(31) 


Ui 


2 

1 


1/15  {  p  -  hu)  +  1/4  pbj^ 


(32) 


and 


I.  =  p  -  b  V 


It  can  be  shown  that 


is  positive. 


(33) 


Again,  Eq.  (31)  is  the  equation  of  motion  for  the  radial  coordinate 
of  a  fictitious  harmonic  oscillator  in  two  dimensions,  but  with  the  spring 
constant  varying  with  time.  If  ^  and  tj  are  Cartesian  coordinates  for 
this  fictitious  oscillator,  they  satisfy  equations  of  motion  of  the  form 


r  ^ 

- 


at  ^  ^  t 

e  +  Wj  ]  ^ 


(34) 


The  substitution 

y  ~ 

I 

converts  Eq.  (34)  to 
d 


_L  jd_ 

y  dy 


y  "cfy 


at 


4  +  [  4 


+  4 


“1 

~z 

a 


y 


=  0 


(35) 


(36) 


« 
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This  equation  has  the  solutions 


4 


2  cj  ^ 


+  ik 


y) 


(37) 


where 


k  = 


(38) 


and  is  the  Bessel  function  of  inriaginary  order 

2^0  _  "  u,„y  2n+ik 

■^ik  n!  '(n+lk)  ! -  '“5“^ 

n=0 


(39) 


The  fictitious  coordinates  ^  and  r\  are  each  a  linear  combination  of 
and  J  and  0  is  then  given  by 


0 


=  v^?T7 


(40) 


That  the  Bessel  functions  enter  the  solution  of  the  small  angle 
problem  is  well  known.  In  fact,  the  principle  difference  between  these 
results  and  those  of  reference  1  is  that  the  frequency  is  different. 


For  early  enough  times,  the  variable  y  is  small,  and  only  the 
first  term  in  the  series  (39)  need  be  considered.  Then,  apart  from 
constant  factors,  the  solutions  (37)  are 


+  ik 

=  y- 


2  wj 
a 


+  i  Wj  t 

e 


(41) 


Of  course,  the  early  time  solution  is  only  valid  provided  the  initial 
"gyroscopic"  motion  is  confined  to  small  angles,  as  was  discussed  in 
the  previous  section.  For  this  case,  =  hp/Z,  and  Eq.  (41)  is  correct. 

It  is  also  worth  noting  at  what  time  the  aerodynamic  forces  begin 
seriously  to  change  the  motion  from  its  initial  "gyroscopic"  form.  Taking 
the  first  two  terms  of  the  series  (39)  leads  to,  apart  from  a  constant 
factor 


t  r,  1  /“a  ,2  - 

^  ®  ^  -  TTIk  <—  )  ^ 


(42) 


where  we  have  used  the  fact  that 


(o  y  =  cj. 


(43) 


Of  course,  u  depends  on  the  time.  We  now  examine  two  cases. 

3i 


Case  1,  «  a  ,  or  rotation  rate  small  compared  to  density  rate. 

In  this  case  k  is  small  compared  to  unity,  and  the  amplitude  of  the  motion 
is  beginning  to  be  seriously  changed  when 

cj  «  Q  (  w,  «  a  )  (44) 

3l  X 


Case  2.  »  a,  or  rotation  rate  large  compared  to  density  rate. 

In  this  case,  k  is  large,  and  the  phase  of  the  motion  begins  to  be  seriously 
changed  when 

2  2 

«  a  “  k  «  2  Wj  a  (phase,  Wj»  a)  (45) 

However,  the  amplitude  of  the  motion  is  not  yet  seriously  changed  at  this 
point.  This  results  from  the  fact  that,  for  large  k  ,  the  terms  in  the 
series  (39)  up  to  n  «  k  approximate  the  series  for  the  exponential 


If  the  n  is  dropped  in  the  factors  (n  +  ik),  one  gets  exactly  this  result, 
which  has  no  change  in  amplitude.  The  correction  is  of  order 


i? 


9 


and  the  amplitude  of  the  motion  is  seriously  changed  when 
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w  ak  «  2cjj  (amplitude,  »  a)  (46) 


We  may  combine  the  results  (44)  and  (46)  into  the  statement  that  the 
amplitude  is  seriously  changed  when  the  air  density  is  high  enough  to  make 

u)  ^  +  4  (amplitude)  (47) 

a 


Finally,  let  us  examine  the  motion  at  late  times,  when  lo  is  large 
compared  to  uj  and  a.  The  asymptotic  form  of  the  Bessel  function  leads 
to 


J. 


ik 


[  1- 


1 

nr 


a 


2 

a 


] 


2^1  2 
“l  “ 


u>  a 
a 


(48) 

- 


Remember  that  increases  exponentially  with  time, 


dw 


dt 


U) 

a 


(49) 


The  "instantaneous  frequency"  of  the  motion  is  the  rate  at  which  the 
argument  of  the  cosine  in  Eq.  (48)  increases,  or 

2^1  2 
1  ^  "T?)  ° 

"instant,  w"  =  u  +  -  (50) 

aw  CO 

a 

Except  for  the  term  in  q^,  this  frequency  is  the  same  as  that  given  by 
Eq.  (22),  provided  the  radical  there  is  expanded  for  large  (Remember 

the  definition  (32)).  We  may  say  that  the  term  is  an  effect  of  the 
changing  "spring  constant". 
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We  also  see  from  Eq.  (48)  that  the  amplitude  of  the  oscillations 
falls  off  as 


amplitude  ~ 


(p  .  )  4 


(51) 


In  fact,  by  comparing  the  initial  and  asymptotic  forms  of  the  Bessel 
function,  one  can  find  the  ratio  of  the  amplitude  to  the  initial  amplitude 


amplitude 
initial  amplitude 


TT  W 


cosh  ( 


(ik) 


2  7ru> 


(■ 


1/2 

1 

(52) 


(See  Jahnke  and  Emde,  formula  f,  p.  11,  for  the  factorial). 

Again,  it  should  be  borne  in  mind  that  this  ratio  is  accurate  only 
if  the  initial  gyroscopic  motion  is  confined  to  small  angles. 


V.  DISCUSSION  OF  THE  LARGE  ANGLE  CASE  BY  USE  OF  THE  ADIA- 
BATlC  invariant. 


The  results  of  Sections  III  and  IV  depend  on  the  assumption  that  the 
angle  9  is  small.  This  angle  does  become  small  eventually,  deep  in 
the  atmosphere,  but  it  can  be  large  initially.  In  this  section  we  shall  use 
the  adiabatic  invariant  to  discuss  the  large  angle  case.  This  treatment 
will  depend  on  the  assumption  that  the  initial  rotational  rate  is  large 
compared  to  the  density  rate  a. 

We  return  to  the  exact  equations,  Eq.  (18)  or  (19).  Initially, 
outside  the  atmosphere,  the  torque  is  negligible.  In  this  case,  Eq.  (18) 
can  be  integrated  once,  giving 

,  (3  -  b  V  cos  0)^  *  4.  -  V  /C2V 

0  +  — - n - =  constant  =  h  (53) 

sin  9 
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2  2 

The  constant  h  is  the  initial  value  of  p  +  q  {see  Eq.  (9))>  the  square  of 
the  angular  velocity  about  the  two  axes  in  the  plane  perpendicular  to  the 
symmetry  axis.  h  could  be  called  the  "tumble"  angular  velocity, 

while  u  is  the  spin  angular  velocity.  h  can  have  any  non-negative 
value,  independent  of  v. 

Eq.  (53)  can  be  simplified  by  the  substitution 

cos  0  =  u,  (54) 


which  leads  to 
2 


ii  =  h  (1-u^)  +  (P  -  bt'u)^ 


This  is  the  energy  equation  for  a  particle  moving  in  a  quadratic  potential, 
and  therefore  has  harmonic  solutions.  The  frequency  S2  q  of  the  oscil¬ 
lations  in  u,  and  0  ,  is  given  by  the  coefficient  of  the  quadratic  term  in 
u,  and  is 

^  +  h  (56) 


Fig.  3.  Relation  of  Various 
Parameters. 


Fig.  3  is  a  vector  diagram  of  the  angular  momenta,  divided  by 
the  moment  of  inertia  A,  for  the  present  case  of  negligible  torque. 

The  total  angular  momentum  l  constant  in  magnitude  and  direction. 
The  constant  p  is  the  projection  of  l/a  onto  the  Z-axis.  bt*  is  the 
angular  momentum,  divided  by  A,  about  the  symmetry  axis,  and  Vh  is 
the  component  of  angular  momentum,  divided  by  A,  perpendicular  to 
the  symmetry  axis.  The  apex  a  of  the  right  triangle  formed  by  a  » 
bv,  and  h  revolves  around  tHe  circle  C  with  angular  velocity 

«o- 


Now  let  the  missile  begin  to  enter  the  atmosphere.  The  torque 
term  in  Eq.  (19)  builds  up  exponentially  like  e .  But  if 


»  Q, 

o 


(57) 


the  change  in  the  force  function  in  Eq.  (19)  per  cycle  2  tt / S2^  will  be 
small,  and  the  action  1, 


.2 

I  =  f  e  d  e  =  f  e  dt,  (58) 

cycle  cycle 


will  be  approximately  constant.  We  shall  evaluate  this  integral  shortly. 

First,  however,  note  that  at  some  altitude,  which  we  shall  call  the 
critical  altitude,  the  torque  term  in  Eq.  (19)  will  become  of  about  equal 
importance  to  the  rotational  term.  At  the  critical  altitude,  the  aero¬ 
dynamic  frequency  is  about  equal  to  S2  ^/Z,  and  the  motion  has  just 
become  seriously  changed  by  the  aerodynamic  forces.  Above  the  critical 
altitude,  the  frequency  of  the  0  variations  is  close  to  fl  q,  while  below 
the  critical  altitude  the  frequency  increases,  following  and  approaching 
2  wa.  Thus,  if  the  adiabatic  condition  is  satisfied  initially,  Eq.  (57),  it 
will  be  satisfied  at  all  times,  since  the  density  rate  a  is  (approximately) 
constant.  The  action  will  therefore  be  approximately  constant  at  all 
times,  both  above  and  below  the  critical  altitude.  Using  this  fact,  we 
can  obtain  a  relation  betw'een  the  initial  rotation  rates,  above  the  critical 
altitude,  and  the  amplitude  of  the  small  oscillations  at  altitudes  somewhat 
lower  than  the  critical  altitude.  To  do  this  we  have  to  evaluate  the  action 
integral  for  the  initial  gyroscopic  motion  and  for  the  asymptotic  small 
amplitude  motion. 

From  Eq.  (53),  the  initial  action  is 


I.  =  /  e  dO  =  2  /  de  \/ a  -  ~  (591 

“  cycle  Sj  e 


where  0  ^  and  ©2  are  the  two  zeroes  of  the  radical.  The  integral  can 
be  simplified  by  changing  the  variable  to  u  =  cos  0  and  can  be  evaluated 
with  the  help  of  formula  (380.  311)  of  Dwight,  p.  71.  The  result  is 
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hr,  =  .  lb  M  )  1 


(60) 


where  |P|  ,  |bi/|  )  means  the  larger  of  |p|  and  |bi/|  .  (This 

choice  corresponds  to  whether  the  tip  of  the  missile,  in  its  gyroscopic 
motion,  does  or  does  not  encircle  the  Z-axis.  ) 

At  altitudes  somewhat  lower  than  the  critical  altitude,  the  aero¬ 
dynamic  forces  dominate  the  motion,  forcing  it  to  small  angles.  The 
form  of  motion  is  then  given  approximately  by  the  results  of  Sec.  Ill, 
in  particular  by  Eq.  's  (24)  and  (25).  Thus  the  final  action  is  approxi¬ 
mately 


,  2 

E.  =  /  0  dt 

fin  , 

cycle 


+  7r/4  u) 

2  / 

t=-jr/4  u) 


4  p  ^  cos^  (2u)t) 
X  +  p  sin  (2  wt) 


dt 


(61) 


Again,  the  substitution  u  =  sin  (2wt)  and  the  same  formulae  in  Dwight, 
enable  one  to  evaluate  this  integral,  with  the  result 


Using  Eq.  (25)  and  the  fact  that  -  X,  we  rewrite  this  equation  as 


E.  =  TT  w  0 

fin 


TT  1  P  -b  I' 


(62) 


Equating  and  we  obtain  an  equation  expressing  the  mean- square 

angle  of  attack  (at  altitudes  somewhat  less  than  the  critical  altitude)  in 
terms  of  the  initial  rotation  rates 

^  =  if  +  2  -  ^(  |p|Jbi;|)J  (63) 

For  to  here,  one  may  use  Eq.  (22),  or  simply  replace  to  by  w^,  since 
the  corrections  are  small  at  later  times.  Eq.  (63)  shows  that  the  mean- 
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square  angle  is  larger  for  larger  initial  rotation  rates. 

We  may  ask  whether  Eq.  (63)  agrees  with  the  result  of  Sec.  IV  in 
the  case  where  the  initial  gyroscopic  motion  is  confined  to  small  angles. 
In  this  case,  (3,  bv,  and  q  are  all  approximately  equal,  differing  only 
by  about  Qg  hv/2,  where  6q  is  representative  of  the  initial  small  angles. 
Equation  (63)  then  becomes,  within  a  factor  of  two  or  so 

?  =  _kii  (64) 

O  CO  '  ' 

a 


This  is  to  be  compared  with  Eq.  (52),  evaluated  for  the  case  where 
<oi/a  is  large,  and  cjj  =  h  v /Z  (see  Eq.  (32)  )  since  p  «  b  This 
evaluation  leads  to  the  same  result  (64). 


VI.  THE  CASE  OF  LARGE  INITIAL  ANGLES  WITH  SLOW  ROTATION 

The  only  case  not  covered  in  previous  sections  is  that  where  the 
initial  rotation  rates  are  small  compared  to  the  density  rate  a  but 
initial  angles  are  large.  We  shall  now  give  an  approximate  discussion  of 
this  case. 

If  the  rotation  rates  are  small  compared  to  a,  they  may  be  ignored 
entirely,  except  insofar  as  they  determine  the  angle  of  attack  the  missile 
will  have  when  it  arrives  at  the  altitude  where  the  aerodynamic  forces 
first  become  appreciable.  Thus  we  may  consider  Eq.  (19)  without  the 
rotational  term.  Let  us  consider  the  torque  function  T(0  ).  It  will  have 
a  shape  something  like  that  indicated  in  Fig.  4.  Up  to  an  angle  0  7r/2, 

a  linear  fit  gives  a  fairly  decent  approximation.  Therefore,  if  the  initial 
angle  is  not  greater  than  about  7r/2,  we  may  use  the  "small-angle”  result 
(52)  evaluated  for  the  case  where  wj/q  is  small,  namely 


amplitude 
initial  angle 


TT  CO 


(65) 


If  the  initial  angle  is  larger  than  about  tt/Z,  the  amplitude  (at  the 
lower  altitudes)  will  be  larger  than  this  formula  implies,  because  the 
first  swing  round  of  the  missile  will  take  longer  than  it  would  if  the  torque 
continued  to  follow  the  linear  fit.  In  fact,  if  the  initial  angle  is  almost 
TT,  it  will  take  a  very  long  time  for  the  first  swing  round.  However,  such 
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cases  are  fairly  improbable. 


Fig.  4 


To  get  further  information,  one  can  resort  to  numerical  solution 
of  Eq.  (19). 
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